The evolution of adult pollen feeding did not alter postembryonic growth in Heliconius butterflies

Abstract For many animals, the availability and provision of dietary resources can vary markedly between juvenile and adult stages, often leading to a temporal separation of nutrient acquisition and use. Juvenile developmental programs are likely limited by the energetic demands of many adult tissues and processes with early developmental origins. Enhanced dietary quality in the adult stage may, therefore, alter selection on life history and growth patterns in juvenile stages. Heliconius are unique among butterflies in actively collecting and digesting pollen grains, which provide an adult source of essential amino acids. The origin of pollen feeding has therefore previously been hypothesized to lift constraints on larval growth rates, allowing Heliconius to spend less time as larvae when they are most vulnerable to predation. By measuring larval and pupal life‐history traits across three pollen‐feeding and three nonpollen‐feeding Heliconiini, we provide the first test of this hypothesis. Although we detect significant interspecific variation in larval and pupal development, we do not find any consistent shift associated with pollen feeding. We discuss how this result may fit with patterns of nitrogen allocation, the benefits of nitrogenous stores, and developmental limitations on growth. Our results provide a framework for studies aiming to link innovations in adult Heliconius to altered selection regimes and developmental programs in early life stages.


| INTRODUC TI ON
Life-history theory predicts that resource partitioning among life stages is optimized by both intrinsic life-history trade-offs and extrinsic environmental effects on survival and reproduction (Partridge et al., 2005;Reznick, 2010;Roff, 1992;Stearns, 1992). In holometabolous insects, ecological and biological differences between developmental periods can cause resource intake and utilization to become separated by life-history transitions, making mechanisms of resource budgeting particularly important (Ahlström, 2011;Boggs, 1981;Istock, 1967). As the main foraging stage, larvae must therefore sustain immediate demands for growth and basal metabolism while building the reserves necessary for the pupal and adult stages (Boggs, 2009;Hahn, 2005).
Lepidoptera provide excellent opportunities to study how variation in resource partitioning can shape life-history traits because they exhibit a variety of foraging ecologies as both larvae and adults (Slansky & Scriber, 1985;Swanson et al., 2016). As a resource, dietary nitrogen plays a crucial role in limiting developmental growth (Mattson, 1980) and constraining reproductive output (Fischer et al., 2004;Cahenzli & Erhardt, 2013;Swanson et al., 2016;Espeset et al., 2019). Since the majority of butterflies and moths feed on nectar (Slansky & Scriber, 1985), which is a poor source of amino acids (Baker & Baker, 1986), the nitrogenous requirements of adult butterflies and moths are largely provided for by nitrogen gathered during the larval phase from host-plant tissue (O'Brien et al., 2002). One notable exception are Heliconius butterflies, which actively collect and digest pollen as a consistent source of amino acids during the adult stage (Gilbert, 1972;Young & Montgomery, 2020). The ability to pollen-feed is a derived trait in Heliconius, marked by distinctive flower-handling behaviors Penz & Krenn, 2000), foraging strategies (Gilbert, 1975), and mouthpart modifications (Krenn & Penz, 1998) absent in all other related genera. Pollen feeding in Heliconius brings about a substantial increase in expected adult nitrogen intake (Boggs, 1981), creating an evolutionary opportunity for change in patterns of juvenile resource allocation (Boggs, 2009). Gilbert (1972Gilbert ( , 1991 and Boggs (1981) therefore predicted that the evolution of adult pollen feeding might have lifted constraints placed on the growth rates of Heliconius larvae by nitrogen nutritional targets, and that, relieved of these constraints, Heliconius may complete their larval development faster and thus minimize juvenile mortality.
Indeed, previous studies in Heliconiini butterflies have highlighted differences in how pollen-feeding Heliconius and nonpollen-feeding outgroups allocate resources to reproduction during metamorphosis (Boggs, 1981;Dunlap-Pianka et al., 1977), indicating that nitrogen targets in Heliconius larvae may indeed be reduced.
Here, we test one of the main predictions about the opportunity provided by the evolution of adult pollen feeding-namely, that Heliconius larvae will develop faster than closely related species with an exclusively nectar-based adult diet. Specifically, we use a comparative approach to study the duration of larval and pupal development, larval growth curves, and pupal and adult weights of three pollen-feeding Heliconius species and three nonpollen-feeding Dryas iulia, Agraulis vanillae, and Dryadula phaetusa (exclusive nectarivores). Adults of each study species were collected within a 2 km radius of Gamboa and placed in 2 × 2 × 2 m cages to mate.
All stock cages received fresh Lantana flowers daily, and artificial feeders with a 35% sugar solution. Heliconius stocks also received fresh Psiguria and Gurania flowers, and crushed bee pollen dissolved in the artificial feeders. Eggs were removed from host plants (Table   S1) by hand each day, and stored without food until emergence.
First, instar hatchlings were transferred to host plant shoots inside individually labeled cups and fed on their preferred host plants; Passiflora biflora (H. erato, D. iulia, and D. phaetusa), or P. platyloba (H. hecale, H. melpomene and A. vanillae), so as to favor optimal growth rates. In addition, to confirm our growth data were consistent with natural populations, body weight data were compared to data from wild-caught adults collected in 2012 and 2013, measured as de-

| Life history characterization
Larval hatchlings were weighed in groups of five, prior to feeding, with individual weight taken as one-fifth of the combined weight.
Larvae were then weighed individually once every two days to minimize handling-related mortality. Individuals that successfully pupated were allowed to dry for 24 h and were then weighed.
The duration of the larval and pupal periods was recorded in days.

All individual weights were measured on a Sartorius H110 Handy
Analytical Balance, with 0.1 mg resolution. Length of pupal period, but not weight, was also recorded for an additional 125 larvae, reared for other experiments.

| Statistical analyses
All analyses were performed using R v.4.0.0 (R Core Team, 2020) and RStudio v.1.3 (RStudio Team, 2020). Unless stated otherwise, linear and generalized linear models (GLM) were generated with lme4 v.1.1-23 (Bates et al., 2015). To test if larval development is shortened in pollen-feeding species, we first assessed interspecific variation in the duration of the larval period (in days) by treating the data as counts and using GLMs with Quasi-Poisson distributions to test the significance of species. Heliconius versus non-Heliconius comparisons were then made to assess whether observed interspecific variation was primarily due to group differences, by means of mixed effect Conway-Maxwell Poisson GLMs with species as a random factor, using the R package glmmTMB v.1.0.1 (Brooks et al., 2017).
Additionally, linear models were built to test the robustness of the results, reported in Tables S2-S7. The same approach was followed to investigate variation in the duration of the pupal period.
We compared larval survival rates across species and adult foraging habits using the R package survival v. 3.2-13 (Therneau & Grambsch, 2000). Larval body mass measurements were used to reconstruct growth curves and ascertain whether patterns of larval growth varied across species and groups. To better detect differences in growth patterns, we normalized developmental age across species for each individual to a range between 0 (hatching day) and 1 (day before pupation). Body mass was also normalized for each individual to a range between 0 and 1, with 1 equating to the last recorded larval weight. Since not all larvae were weighed on the last day before pupation, the difference in days between their final weight measurement and their final day as larvae was introduced as a random effect in the models. Linear models and GLMs with Gamma distribution and log-link functions were built to compare growth curves, with individual and species as additional random effects.
Models of normalized data with the lowest AIC scores are reported in the main text. Additionally, pairwise comparisons were run using the R package statmod v.1.4.34 (Giner & Smyth, 2016). Full results from normalized growth curve models, as well as equivalent models built on the original raw growth data, are presented in Tables S2-S7. Finally, we detected interspecific variation in pupal and adult body weight using linear models, while linear mixed models with species as a random effect were used for Heliconius versus non-Heliconius comparisons. Linear regressions with species as a random effect were then built to test whether the duration of the larval period correlated with pupal and adult weight.

| Larval development
We tracked the development of 125 larvae, of which 58 reached the 5th and final instar. Survival rates did not vary significantly across species (X 2 5,125 = 8.90, p = .100) or between pollen feeders and nonpollen feeders (X 2 1,125 = 1.00, p = .300). While we found significant variation in the duration of the larval period between species (F 5,52 = 10.30, p < .0001), pollen feeders did not spend less time as larvae than nonpollen feeders (X 2 1,58 = 0.78, p = .379; Figure 1c; Table 1). Instead, interspecific variation in the duration of larval period was primarily driven by D. phaetusa, which had the longest developmental time (Table S2). When normalized to control for variation in the length of the larval period, larval growth curves ( Figure 1a), which capture variation in the dynamics of the larval growth phase, also differ significantly between species (Figure 1b; F 5,57 = 27.15, p < .0001). Although pairwise contrasts were inconsistent across models, growth curves of pollen-feeding Heliconius consistently do not differ from their nonpollen-feeding relatives (F 1,57 = 0.05, p = .829; Table S3).

| Pupal development
Fifty-seven of the tracked larvae pupated successfully, and 46 completed pupal development. The weight of dried, fresh pupae varied across species (F 5,51 = 17.74, p < .0001), but this variation was not explained by presence or absence of adult pollen feeding (F 1,3.98 = 0.04, p = .845; Table S4). Heliconiini pupal development spanned 6 to 9 days, with A. vanillae, H. erato, and H. melpomene pupae developing significantly faster than H. hecale, and significantly slower than D.
iulia (F 5,119 = 17.51, p < .0001; Table 2; Table S5). Interspecific variation in the duration of the pupal period was not explained by pollen feeding (X 2 1,125 = 0.38, p = .540).  ). To confirm that our rearing conditions did not bias patterns of growth, we compared our adult weights to those of wildcaught butterflies from the same population. We found that wild butterflies were lighter than caged ones (F 1,92 = 20.65, p < .0001), but the ranks between species were conserved (r(4) = 1.00, p = .003).

| DISCUSS ION
By evolving to utilize a novel source of protein in the adult stage, Heliconius butterflies are hypothesized to have had an evolutionary opportunity to shorten the time spent as vulnerable larvae (Gilbert, 1972(Gilbert, , 1991. Collecting pollen as adults, which provides a rich source Heliconius. First, pollen-feeding H. charithonia and H. cydno eclose with less total abdominal nitrogen, a proxy for larval reproductive investment, than the nonpollen-feeding Dryas iulia (Boggs, 1981).
Second, the ovaries of H. charithonia females are also smaller and contain fewer total oocytes than those of D. iulia females (Dunlap-Pianka et al., 1977). These data strongly suggest that pollen feeders earmark relatively fewer larval reserves to reproduction, and if this energetic requirement were reallocated to the adult stage, this could favor the evolution of faster larval development in Heliconius. Below, we discuss three plausible alternative hypotheses that may explain why the energy diverted away from reproduction in Heliconius is not reallocated to larval growth rates.
First, Heliconius butterflies may use the larval reserves not invested in reproduction to strengthen their chemical defence (Cardoso & Gilbert, 2013). If so, the total nitrogen nutritional targets may remain the same for all Heliconiini larvae, and therefore limit growth rates across the board. Heliconiines are chemically F I G U R E 1 Larval growth curves and duration of larval development in pollen-(red tones) and nonpollen feeders (blue tones). (a) Variation in growth patterns and duration of larval period in raw larval weight data for 58 fully tracked individuals. Trend curves are generated with ggplot2 loess function and have no statistical purpose. (b) Variation in growth patterns and duration of larval period in normalized larval weight data for 58 fully tracked individuals as used in statistical comparisons. Growth patterns differ significantly between species (F 5,415 = 3.202, p = .008) but not by adult foraging strategy (F 1,417 = 0.292, p = .589). (c) Duration in days of the larval period per species. Duration of the larval period did not differ between pollen-and nonpollen feeders (X 2 1,58 = 0.398, p = .528). Asterisks denote significant interspecific contrasts at p < .0001 in posthoc comparisons. (d) Phylogenetic relationships between study species, adapted from Kozak et al. (2015) defended by cyanogenic glycosides, both sequestered from their Passifloraceae host plants, and biosynthesized de novo from amino acids (Davis & Nahrstedt, 1987;Engler-Chaouat & Gilbert, 2007 -Pianka et al., 1977) and are affected by the quality of larval diet (Dunlap-Pianka, 1979). Similarly, for at least some Heliconius species, cyanogenesis is unaffected by pollen deprivation during the first 20 days after emergence, but declines thereafter (Cardoso & Gilbert, 2013), suggesting that cyanogenic metabolism also exploits larval reserves during early adult life. These data indicate that larval-derived resources in Heliconius likely have a large role in supporting both reproduction and toxicity during the first weeks of adult life. This scenario appears to be at odds with cited evidence of comparatively less abdominal nitrogen and reproductive tissue in Heliconius than in nonpollen-feeding Dryas iulia (Boggs, 1981;Dunlap-Pianka et al., 1977). However, organism-wide nitrogen contents in Boggs (1981, personal communication), are comparable between H. cydno and D. iulia (but not H. charithonia), suggesting that nitrogen not allocated to reproduction can be redirected elsewhere following adaptive pollen use, although this may not always be the case. More research is needed to track larval-and adult-derived nitrogen in the products of reproductive and cyanogenic processes throughout the first month of the adult stage.
Third, even if Heliconius larvae do have a smaller nitrogen quota to fill, nitrogen accumulation may not be the limiting factor for faster larval development (Dmitriew, 2011). The larvae in this study were grown in their native environment under ambient conditions, with abundant supply of their preferred host plants and minimum incidence of disease. Therefore, we have likely observed the higher range of each species' optimal growth rate. Our results may simply reflect the developmental constraints that cap all Heliconiini larvae under ideal growth conditions, and a faster development would not be possible without substantial costs (Arendt, 1997). For example, it may be that the developmental tasks shared by all Heliconiini, such as the accumulation of other nutrients, or the time and resources required to form the precursors of adult tissues or synthesize certain metabolites, are more time-limiting than the formation of nitrogen reserves.
In conclusion, the simple prediction that increased quality of adult diets will relax constraints on Heliconius growth rates during vulnerable larval and pupal stages is not supported by our results.
Indeed, current available evidence suggests that pollen feeding may have primarily relieved constraints on adult, rather than larval, development, allowing Heliconius to dramatically extend their lifespan (Ehrlich & Gilbert, 1973;Gilbert, 1972), sustain life-long oogenesis (Dunlap-Pianka et al., 1977), and potentially free nitrogenous resources for use in enhanced chemical defence (Cardoso & Gilbert, 2013;Young & Montgomery, 2020). To fully understand the consequences that the evolution of pollen feeding had on life-history traits in Heliconius, it will be important to investigate trends in re-

ACK N OWLED G M ENTS
The authors thank Oscar Paneso, Cruz Batista Saez, and Wyatt Toure for support at the insectaries, and the EBaB lab (Bristol) and Butterfly Ecology and Evolution Research group (Smithsonian) for advice and feedback. They also thank Tom Finch and Lee Dietterich for statistical advice. Finally, the authors are grateful to the Ministerio del Ambiente, Panama for collection permits and facilitating this work.

CO N FLI C T O F I NTE R E S T
We declare no conflicts of interest.

DATA AVA I L A B I L I T Y S TAT E M E N T
Data are available on the Data Dryad Digital Repository, https://doi. org/10.5061/dryad.3r228 0gjn.